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The possible chiral phase transition in high energy heavy-ion collisions may lead 
- ^ to the formation of a disoriented chiral condensate (DCC). However, the existence 

C _ J , of many uncorrelated small domains in the rather large interaction volume and the 

huge combinatorial background make the experimental search for the DCC signal 
a rather difficult task. Wo present a novel method for studying the formation of a 
DCC in high energy heavy-ion collisions utilizing a discrete wavelet transformation. 
(N ■ We find the wavelet power spectrum for a DCC to exhibit a strong dependence on 

^ , the scale while an equilibrium system and the standard dynamical models such as 

. HIJING have a flat spectrum. 

(N 

■ In heavy-ion collisions the rapid expansion of the collision debris in the 

I longitudinal direction may lead to supercooling of the interaction region, and 

' as a result, domains of the "unconventionally" oriented vacuum configurations 

allowed by the chiral symmetry may be formedlJ Detection of this interest- 
ing phenomenon, the so-called Disoriented Chiral Condensate (DCC)lil, would 
C^hI provide valuable information on the vacuum structure of the strong interac- 

tion and the nature of the chiral phase transition. Preliminary theoretical 
I investigations on non-equilibrium dynamics using the classical linear cr-model 

• have found some evidence for the growth of long wavelength Piiftp modes, 

which may indeed lead to domain structure or cluster formationH'B If there 
are many uncorrelated small domains, the integrated probability distribution 
of the neutral pion fraction / emitted from a disoriented region, predicted to 
be P{f) = l/2y^,a would become Gaussian. This makes the experimental 
search for the DCC signal a rather difhcult task. In order to disentangle the 
DCC domain structure in high energy heavy-ion collisions, we propose a new 
method which emphasizes not only the behavior of the probability distribution 
in the full phase space region but also its fluctuation in rapidity rj or azimuthal 
angle 4>. In other words, one needs to study the "local" properties of the dis- 
tribution in phase space if the DCC clusters are "localized" objects in phase 
space. 

A novel method is a multiircsolution analysis performed by a discrete 
wavelet transformation (DWT)Ij which has been found effective in systemat- 
ically detecting structures of various scales in turbulence, and multiparticle 
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productionsB'El For any one-dimcnsional sample of a point distribution f{x), 
such as the rapidity distribution of the neutral pion fraction in the interval 
[0, 1] with resolution Ax, can be represented as a histogram of 2'^ bins where 
J = mod(|lnAa;|/ln2) + 1 

2'-l 

/(x) = /('^)(a:)= ^ (1) 

where ipjkix) are the piecewise constant functions called the mother functions. 
The histogram function at the next finer resolution is the average of adjacent 
bins at the finest scale. Evidently some detail is lost compared to (|l|). The 
difference between (|l|) and /^'^"^^ can be fully expressed in terms of the diffence 
functions '0j_i_/j(a;) which are stairwell functions called father functions. The 
difference amplitudes or one defines the DWT coefficients 

fjk^2^ J f{x)i;,kix)dx, (2) 

where the basis functions tpjkix) have a unique property: they are related 
to each other by scaling and translation: ipjkix) = ip{2^x — k). The index 
j denotes the resolution scale and k the position at each scale. The idea of 
multiresolution analysis is to find representations of the sample function f{x) 
at various scales. Unlike the Fourier transformation which requires the in- 
formation in whole physical space, the wavelet SSD requires only the local 
information in space. Roughly speaking, the multiresolution analysis is a jour- 
ney down (or up) a hierarchy of scales and to view the "world" in different 
"eyes" with magnifying or reducing glasses. 

The central idea of using the wavelet analysis is to attempt an unbiased 
separation of the DCC physics from the various "backgrounds" which cannot 
be otherwise easily eliminated. The usual correlation method is not so effec- 
tive here because the different scale correlations are mixed together in the same 
phase space. The discrete wavelet transformation (DWT), on the other hand, 
has an extra dimension - the scale which can be used to separate these cor- 
relations. The main source of the "backgrounds" is the fluctuations at small 
correlation scale due to pion radiation from the uncorrelated regions in space- 
time, which can coexist with a DCC domain. The background can also come 
from the partial thermalization of the system leading to thermal correlations. 
We have calculated the wavelet power spectrum for rapidity distribution of the 
neutral pion fraction simulated using the classical linear a modelB One defines 
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Figure 1: The power spectra of the wavelet coefficients for different dynamical scenarios. 



the power spectrum with respect to the wavelet basis at each j as 

2^—1 
k=Q 

Pj is then the power of fluctuations on length scale It is also the mean 

square of the probability distributions of the wavelet coeflicients defined in 
. The result is shown in Fig. 1 . For the random samples and the equilibrium 
evolution of the classical field where there is no DCC production mechanism, 
the wavelet power is fiat indicating there is no sensitivity on the change of the 
scale. For the DCC samples, the power starts to build up as j goes down below 
jd 4 and rises rather quickly, departing distinctively from the equilibrium 
and the random cases. This suggests the existence of the DCC clustering with 
a typical size of Arjd = 27/niax/2-''' — 0.8 units in rapidity. Also plotted is 
the power spectrum from HIJING Monte Carlo data which also features a flat 
distribution. 
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